Dynamics of Polymer Decompression: Expansion, Unfolding and Ejection 
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The dynamics of polymer decompression from a compact state to swollen conformation can be 
formally described as nonlinear diffusion. We discuss two basic examples: (i) the expansion, or 
unfolding from a compact state, and (ii) the ejection of a compressed polymer through a pore. The 
problem can be solved exactly for case (i), but not for case (ii). Even in such situations, a scheme 
called uniform approximation is shown to be useful to get a physical insight involved. Its application 
to the case (ii) is able to account for conflicting numerical data in a consistent way. 
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Macromolecules assume compact conformations in cer- 
tain situations. Examples include DNA in living cells, 
proteins in native states and other polymers in poor 
solvent conditions or under a compression field. Once 
released from the condition, these polymers expand to 
swollen coiled state, which is characterized by developed 
fluctuations. This expansion process is of interest in two 
different contexts. First, this process is relevant to coil- 
globule transitions, thus, regarded as a fundamental topic 
in polymer science [U [2] . Unfortunately, there are com- 
paratively fewer past studies on this expansion process 
than on the folding (coil to globule) process. [3J |H |S] 
The second case of interest is encountered in the field 
of confined polymers [5] . A recent advance in nanoscale 
fabrications and single-chain experiments allows one to 
manipulate and observe individual polymers, thereby fa- 
cilitating a number of potential applications in biological 
as well as nanoscale sciences [7] [5J EE EH EQ H2] • The aim 
of the present Letter is to provide a unified framework to 
describe decompression processes engaged in such situa- 
tions. 

Naively, the problem is analogous to the diffusion of 
molecular gases, which are initially confined in a finite- 
size box. For our case, the essential difference from 
this simple example lies in the connectivity of monomeric 
units into a string, in which a small entropy does not play 
any major role. By taking such polymeric natures cor- 
rectly into account, it is naturally formulated as a nonlin- 
ear diffusion problem, which serves as a basis for various 
decompression processes such as (i) the expansion, i.e., 
unfolding, and (ii) the ejection of a geometrically com- 
pressed polymer from a narrow pore (Fig. [lj . 

In the case of polymer expansion (case (i)), a symme- 
try allows for an analytic solution. However, this is a 
rather special case, and the nonlinearity in the partial 
differential equation usually requires a numerical calcu- 
lation. In such a case, we propose a scheme uniform ap- 
proximation, in which the nonuniformity of the segment 
distribution is totally neglected. It allows one to obtain 
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FIG. 1: Schematics of [top] expansion (unfolding) and [bot- 
tom] ejection processes. 



an approximate solution easily, which can yet be rea- 
sonably compared with the exact solution, provided that 
mechanisms of the driving force and the dissipation are 
correctly identified. We apply this scheme to the case 
of polymer ejection (case (ii)) and demonstrate a very 
good agreement with the numerical simulations reported 
so far. 

Consider a compressed globular chain, the size Rq 
of which is substantially smaller than the relaxed size 
Rp = aNQ in bulk (a, No are the size, the number of 
monomers and v is the Flory exponent). The corre- 
sponding average volume fraction </>o ~ N§o? jR\ <C 1 
is assumed to be small. Here and in what follows, 
we assume that isotropic compression, i.e., in a sphere, 
but generalization to anisotropic compression is straight- 
forward. When given an escaping route, the polymer 
starts to be decompressed toward the relaxed coil. The 
time evolution of volume fraction </>(r) obeys the con- 
tinuity equation; d<fi(r)/dt — — V • j(r) with the flux 
j(r) = — L(r)[V/x(r)]0(r). In order to obtain the chem- 
ical potential /x(r) = df{r)/d<j){r) and the transport co- 
efficient L(r), we utilize the fact that the compressed 
polymer resembles a semidilute solution, which amounts 
to say that the compressed polymer is viewed as a dense 
piling of blobs of size £(r) ~ a0(r)"/( 1— 3 "'. Such a situ- 
ation is called a strong confinement (SC) in Refs. [131114] 
(see also eq. §5§). Then, the excess free energy density 
due to the compression and the dissipation mechanism 
are, respectively, evaluated by assigning ~ fcgT (thermal 
energy) and the Stokes friction per blob [15]: f(r) ~ 
fcsT/£(r) 3 and L(r) ~ g( r )/( T l£,( r )) where n is the solu- 
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tion viscosity and g is the number of monomers per blob 
(f ~ aflr"). This leads to j = -a 2 /T <Kr) i//(3l/ ~ 1) V^(r) 
with the monomer scale time constant To — rja 3 /{kgT). 
The same expression is obtained from slightly different 
route, which utilizes the semidilute solution analogy in 
a more transparent way: j = — D c {r) V0(r), where 
D c (r) ~ ksT / (r]£(r)) is the cooperative diffusion coeffi- 
cient [15]. Using this flux and the continuity equation, 
we arrive at the following: 



9 t 0=(a 2 /r o )V[(0(r)) m V(^(r))], 



(1) 



with m = v/(Zv — 1) |28j . This equation constitutes a 
basis for the decompression process and should be solved 
with appropriate initial and boundary conditions speci- 
fied for each case of interest. 

a. Expansion and unfolding There are several ways 
to achieve decompression. The simplest is to switch off 
the external compression field abruptly at t = 0. Then, 
the polymer will expand isotropically toward the relaxed 
coil(Fig. IT]). In the context of the unfolding by the tem- 
perature (more generally, solvent quality) jump across 



Tfi n ), the initial compact 



the 9 temperature (T i? 
state can be viewed as a dense piling of thermal blobs 

1 /2 

of size £th = a/Kim| = ag t h , where <; = (T — 6)J9 is 
the reduced temperature. Hence, the initial size of the 

polymer is R ~ Cth(N /9th) 1/3 = a< /3 /|Q m | 1/3 . 

Let us first make eq. (JlJ dimensionless by introduc- 
ing a scaled coordinate (radial component) f = r/Rp, 
volume fraction <p(f) = (f>(r)/(f>F and time t = t/r, 
where 4>f = o, 3 Nq/R f denotes the monomer volume 
fraction in a relaxed state. By setting the time unit 
t = Tex — tqNq 1 ', we can rewrite eq. |lj as 



(2) 



The characteristic time for expansion r ex depends only 
on the relaxed coil size Ftp and not on the initial size. 
Therefore, it has the same scaling form as the longest re- 
laxation time in the relaxed state. The solution reads 125 



r,t) = c R(i)-' 6 \-{f/R{t)} 2 



1/m 



(3) 



where c = T[5/2 + l/m\/(ir 3 / 2 T[l/m + 1]) with the 
Gamma function, T, is the normalization factor and R 
corresponds to the overall chain size, which evolves ac- 
cording to 



R(t) = Ro 1 + t/f, 



(4) 



l)/(9v-2) and the time 



' ?1/Q 7(2c re ). Plots of the solu- 



with the exponent ctQ = (3v 
constant f/°] = aomR^ 

tion are shown in Fig(2]3 along with the results of nu- 
merical calculation carried out using the initial condi- 
tion of a step-like distribution. In the strict sense, so- 
lution ^ is valid only for the initial condition <fio — 
c -Rg~ 3 [l — (r/Ra) 2 ] 1 /™. Nevertheless, the figure shows 



that the dependence on the initial distribution is not 
significant. The analytical results correspond well with 
those of the numerical calculation, particularly for the 
mean size of a chain R m = J rij>(r)dr/ J <f>(r)dr (~ R) 
(Fig§i). 
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FIG. 2: Solutions of nonlinear diffusion equation (solid lines) 
(with v = 3/5) together with numerical results for initial con- 
dition of step-like distributions (points). 

b. Uniform Approximation Except for fortunate 
cases as above, one has to rely on the numerical compu- 
tation. In such cases, it would be desirable to find out an 
approximation scheme allowing for a physical insight into 
the process. In this spirit, we propose uniform approx- 
imation, in which the spatial distribution of monomers 
are regarded as uniform. Note that almost the same at- 
titude was adopted to study the stretching dynamics of 
a tethered polymer (monoblock picture) |22j . 

Let us illustrate it by applying to expansion dynam- 
ics. In the expansion process, monomers are assumed 
to be uniformly confined in a spherical cavity of size 
R(t), which is the size of the polymer chain at the mo- 
ment(Fig. [TJ. The volume fraction <fi(t) ~ a 3 N /R(t) 3 is 
connected with the blob size f(t) ~ a</>(t) l//(1_3l/) . Then, 
the excess free energy of the confinement is written as 



AF/(k B T) ~ R(t) 3 /atf * [a/R(t)] 
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A nonlinear dependence on the chain length: A_F ~ iv 
for a chain with v = 3/5 (in good solvent) and AF ~ Nq 
with v = 1/2 (in 8 solvent) is a characteristic of the SC 
regime [T3J H3[ . By noting the velocity gradient Vw ~ 
the overall dissipation is evaluated as 



TS(t) ~ r)R(t)R(tY 



(6) 



Balancing the rate of the free energy change and the dis- 
sipation, we obtain the following solution 



AF(t) = -TS(t) & R(t) = R (1 + t/ Tmi f 



(7) 



with the exponent a = (3v — l)/(9v) and a time scale 

T ini = T o(Ro/ a NQ^ 3 ) 1 / a - Again, this can be rescaled with 
the characteristic time T ex so that we obtain the same 
form as in eq. Q; R(t) = R(t)/R F = R Q (l + t/f ml ) a 
with the rescaled time fj n j = Ti n i/r ex — (R$) 1 / a . This 
time fj„j designates an initial transient period, which de- 
pends on the degree of the initial compression C com = 
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4>q/4>f as Ti n i = Ccom ■ At t > f im , a sub-diffusion 
like scaling regime R(t) ~ t a follows. All these results 
agree reasonably well with the exact solution except for 
a slight difference in the dynamic exponent. In partic- 
ular, the uniform approximation is able to predict the 
correct characteristic time T ex of expansion. This success 
traces back to the correct evaluation of the driving and 
dissipation mechanisms based on the fundamental length 
scale £. 

c. Ejection Consider a long polymer initially con- 
fined in a spherical cavity of size Rq. Through a small 
hole (with a size comparable to that of monomers), the 
polymer escapes (Fig. [I]). This process of polymer ejec- 
tion, i.e., a polymeric version of effusion, has recently 
attracted considerable attention inspired by the biolog- 
ical problem such as a DNA packed into and unpacked 
out of a phage capsid. From a fundamental point of view, 
there have also been many theoretical and numerical at- 
tempts to understand how a flexible chain is ejected from 
a pore [THl Q~71 HH HH] • Even in this simplest setup, how- 
ever, the situation is yet controversial. It was reported 
that the Monte Carlo estimate of the ejection time r ej 
can be fit well by the formula 
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which may be obtained by assuming AF ~ AFideal ~ 
No [IZ]. One of the latest works criticized eq. ^ by 
emphasizing the importance of the notion of SC, and 
proposed the following formula based on the systematic 
Monte Carlo simulation 118 : 
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iV 1+ 7A M ~ N 1+ »ct> 



l+u ,1/(1-31/) 



(9) 



where Afj, — AF/No is the free energy of the confine- 
ment (eq. ([5])) per monomer, i.e., the chemical potential 
difference as a driving force. However, the origin of this 
scaling is unclear. Indeed, it was proposed as a translo- 
cation time of a polymer driven by strong and constant 
Aji |20] , which would be asymptotically valid in the limit 
of the strong driving force [5T] . In the present case of the 
ejection, A^i is far below such an asymptote and grad- 
ually decreases with the process advanced. Moreover, 
eq. (|9| shows an unexpected opposite trend of the chain- 
length dependence for polymers in good and 6 solvent: 
T^f ~ for v = 3/5 and T M „ jv-i/2 for v= i/ 2 . 

In such a situation, it is pertinent to go beyond the scal- 
ing formula and to clarify the underlying physics. 

In principle, one can numerically solve eq. ([!]) with 
the boundary condition ensuring a vanishing flux at the 
wall except for the pore region. Here, we shall ap- 
ply the method of uniform approximation to gain use- 
ful insight. The free energy is essentially the same 
as that in eq. ([5]), but now the dynamical variable is 
the number of monomers N(t) left inside the cavity, 
while the cavity size Ro remains constant: AF/ (&bT) ~ 
(a/ Ro) 3/{3 "- 1} N(t) 3 "/ < - 3 "-- 1 \ At a particular moment 
during the ejection, only the monomer residing at the 



pore feels a resultant osmotic driving force ~ fc,gT/£(t). 
This affects the other monomers yet interior cavity, which 
would be quantified by a distance dependent response 
function. A linear response theory then connects the re- 
sponse to the correlation function of the concentration 
fluctuation. This consideration indicates that the over- 
all dissipation is dominated by the proximity of the pore 
within the range of the correlation length £(t), where 
there exists a velocity gradient of segments of the order 
~ aN(t)/£(t). Therefore, the dissipation term is evalu- 



ated as TS(t) ~ 77 
tion is obtained as 



N(t)a/£(t) £,(t) 3 . The time evolu- 



AF(t) = -TS(t) N(t) = N Q (1 + t/ny 



(10) 



with the exponent /3 = (1 — 3f)/(2(l — v)) and time 

j(l+")/(l-3») ; 



Nq. The dynamic feature 
in eq. ( 10 ) can be well compared with the recent result 



constant t\ ~ tq<j 



of the stochastic rotation dynamic simulation [19 . This 
osmotic driven process ceases at ti that is determined by 
N(t 2 ) ~ g n,2 = ti + t 2 = T <j>, 



-{v+2)/(3v) N {2+v)/(?,v) 











in which g — (Ro/a) 1 / 1 ' is the number of monomers 
corresponding to the relaxed chain size Ro- After T2, 
there is no driving force, and the chain diffusively escapes 
from the cavity, which takes a time td ~ Rq/Dn with 
the diffusion coefficient D^ of the polymer. Thus, we 
arrive at the formula for the total ejection time r ej — 
T2 + TD — Ti.2 + T£), where the second near-equality holds 
in the situation of our interest, i.e., Rf 3> Rq- 

Let us see how this can be compared with previous sim- 
ulations. We note that T e j was so far measured in a model 
system in which the effect of the induced flow of sol- 
vent, i.e., hydrodynamic interactions, was neglected. To 
cope with this situation, we need to modify the dissipa- 
tion term as TS(t) ~ Tja 3 [N(t)] 2 g(t). This modification 
changes the dynamic exponent to j3 = (1 — 3^) /(3(1 — v)) 
and the time constant to t\ = to~ 



To^o 

D No 



3v) 







Nq and n,2 = 
Identifying the diffusion coefficient as 
ksT I (rjaNo), the diffusion time is evaluated as 
td = To4> 2 ^ 3 Nq^ 3 . In the limit of strong compression 
Ro C iJy with large N and Ro, r e j asymptotically ap- 
proaches the simple scaling form ~ r (-Ro/ a ) 3 ^ with no 
dependence on Ao- However, in most practical situations, 
the two contributions n 2 an d td are mixed. Analyzing 
various reported data in terms of our formula strongly in- 
dicates that the result in ref. [T7] (with small Nq and Ro) 
corresponds to the late-stage diffusion dominated regime. 
In contrast, ref. [18] employing larger No and Ro most 
probably monitors the mixing of the driven-ejection and 
diffusion regimes. The asymptotic regime is realized only 
for extremely long chains. Examples of the analysis (with 
v = 3/5) are shown in Fig. [3] in which the dependence 
of r e j on R shows a clear crossover from the diffusion 
dominated regime (r e j ~ 773 ~ Rq, whose scaling form is 
close to eq. (|8]|) to the asymptotic (r ej ~ ri, 2 ~ Rq), be- 
tween which there is a region where eq. r e j ~ Rq^^ 
looks valid. In addition, scaling plots ofreported ejec- 
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FIG. 3: Comparison of predicted ejection time with that ob- 
tained from previous numerical experiments, (a) Double log- 
arithm plot of r e j as a function of Rq for various system sizes, 
(b) Scaling plot of r e j vs. td for small systems, (c) Scaling 
plot of r e j vs. 7-1,2 + td for intermediate systems. Numerical 
data of T e j in (b) and (c) were extracted from [17] and [18] . 
respectively. 

tion times provide more concrete evidence that ref. [1~7] 
and ref. [18], respectively, monitor the ejection processes 
in the small and intermediate size systems in accordance 
with the above scenario. We obtain very good data col- 
lapse without any fitting parameters (Fig. 3b and c) [27] . 
This supports the suitability of our analysis. 

In conclusion, we have proposed a general framework 
for the decompression dynamics of polymers. Aside from 
a basic equation of the nonlinear diffusion type, uniform 
approximation has been shown to be useful. We applied 
the idea and method to the problems of polymer expan- 
sion, unfolding and the ejection. We believe that, for the 
ejection problem in particular, this provides the first-ever 
consistent description with both the driving force and the 
dissipation mechanism clearly identified. 

We would like to end with some remarks and per- 
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